NATIONAL  ADVISORY  COMMITTEE 
FOR  AERONAUTICS 

TECHNICAL  NOTE  2692 

ON  THE  FORM  OF  THE  TURBULENT  SKIN-FRICTION  LAW 

AND  ITS  EXTENSION  TO  COMPRESSIBLE  FLOWS 

By  Coleman  duP.  Donaldson 

Langley  Aeronautical  Laboratory 
Langley  Field,  Va. 


Washington 
May  1952 


DISTRIBUTION  STATEMENT  \ 
Approved  for  Public  Release 
Distribution  Unlimited 


Reproduced  From 
Best  Available  Copy 


je msqpMtzrz  mssmsm  # 

■Q'ty  fhOO~(0-37 72- 


20000731  139 


NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 
TECHNICAL  NOTE  2692 

ON  THE  FORM  OF  THE  TURBULENT  SKIN-FRICTION  LAW 
AND  ITS  EXTENSION  TO  COMPRESSIBLE  FLOWS 
By  Coleman  duP.  Donaldson 

/SUMMARY 


A  derivation  of  the  form  of  the  incompress ible  turbulent  skin- 
friction  law  for  an  insulated  flat  plate  is  made. in  such  a  way  that  it 
may  be  extended  to  compressible  flows.  The  ratio  of  compressible  to 
incompressible  skin  friction  is  obtained,  and  the  results  are  shown  to 
be  in  agreement  with  existing  experimental  results. 


INTRODUCTION 


The  magnitude  of  the  skin-friction  drag  encountered  by  a  flat  plate 
immersed  in  a  fluid  at  Reynolds  numbers  large  enough  to  insure  turbulent 
flow  has  been  one  of  the  basic  problems  of  aerodynamics.  From  the  theo¬ 
retical  approaches  of  Prandtl  and  von  Karman  (for  a  resume,  see  refer¬ 
ence  l)  and  the  experimental  work  of  numerous  investigators,  principally 
Nikuradse  and  Ludwieg  and  Tillmann  (see  references  2  and  3) ?  much  has 
been  learned  of  the  forms  of  the  turbulent  boundary- layer  velocity  pro¬ 
file  and  the  skin  friction  associated  with  these  forms  in  incompressible 
flows,  although  the  exact  mechanisms  involved  are  still  not  completely 
understood. 

Recently,  the  magnitude  of  the  turbulent  skin  friction  on  a  flat 
plate  at  high  Mach  numbers  has  become  of  great  interest,  and  several 
papers  have  been  written  on  this  subject  presenting  both  theoretical 
treatments  of  the  problem  and  the  results  of  skin-friction  measurements 
at  Mach  numbers  between  1.5  and  3.0.  (See  references  4  to  8.) 

The  agreement  between  these  theories  and  the  experimental  data  that 
exist  is,  in  general,  satisfactory.  The  status  of  the  problem/  however, 
is  such  that  a  simple  physical  approach  to  the  extension  of  the  incom¬ 
pressible  skin-friction  laws  to  the  compressible  case  would  seem  desir¬ 
able.  The  purpose  of  this  paper  is  to  present  such  a  simple  physical 
picture. 
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SYMBOIS 


parameter 


[  |"n(r  -  l)  v0 

IL  k2  »08j 

/ 

skin-friction  coefficient 


/  2tw 

Wo 


roughness  height 

constant  relating  mixing  length  to  y 
mixing  length 
constants 
Mach  number 

definitive  ratio  of  total  shear  stress  to  laminar  shear  stress 


Reynolds  number 


Reynolds  number 


Reynolds  number 


absolute  temperature 
velocity  in  x-direction 
velocity  in  y-direction 
distance  along  surface 
distance  normal  to  surface 
ratio  of  specific  heats 


boundary-layer  thickness 
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[i  viscosity 

V  kinematic  viscosity 

p  density 

t  total  shear  stress 

Subscripts : 

L  conditions  at  edge  of  laminar  sublayer 

o  free-stream  conditions 

lam  stress  produced  by  laminar  action  alone 

turb  stress  produced  by  turbulent  action  alone 

w  wall  conditions 


* 

THE  INCOMPRESSIBLE  SKIN-FRICTION  LAW 


If  the  shearing  stress  due  to  turbulence 
velocity  gradient  is  assumed  to  be  given  by 


in  a  region  subject  to  a 


Tturb 


,2  du 
pz  ¥ 


Su 

Sy 


which  is  the  mixing-length  formula  of  Prandtl  (see  reference  1,  p.  130), 
the  equation  of  motion  to  be  satisfied  by  the  boundary  layer  on  a  flat 
plate  is 


Su  +  v  '  du  _  1  St  _  1  _S_  Su 
Sx  Sy  P  Sy  P  Sy  Sy 


2 


Many  investigators  have  established  that  for  most  purposes 


(2) 


b 
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is  an  adequate  representation  of  the  houndary-layer  velocity  profile 
outside  a  very  thin  laminar  sublayer.  For  a  flat  plate  in  the  absence 
of  pressure  gradient  the  value  of  n  is  approximately  7-  It  is  also 
well-known  that  the  value  of  the  mixing  length  in  the  boundary  layer 
near  the  laminar  sublayer  is  given  by 


i  =  ky 


(3) 


where  k  is  approximately  0.4. 

When  the  preceding  facts  are  used,  the  ratio  of  the  total  stress 
to  the  laminar  stress  at  points  in  the  boundary  layer  near  the  wall 
may  be  found  in  the  following  manner: 


Since 


the  insertion  of  l  from  equation  (3)  and  5u/Sy  obtained  from  equa¬ 
tion  (2),  namely, 


1-n 

du  =  uQy  n 
Sy  1 


yields 


+ 


n+1 

u0k2y  n 

1~ 

nv6n 


(4) 


This  equation  is  consistent  with  the  idea  that  the  total  stress  in 
the  boundary  layer  becomes  principally  laminar  as  the  surface  of  the 
plate  is  approached.  The  extent  of  this  laminar  sublayer  may  be  reckoned 
by  computing  the  value  of  y  where  the  ratio  of  the  total  stress  to  the 
laminar  stress  is  given  by  some  definitive  ratio.  When  this  definitive 
ratio  of  total  stress  to  laminar  stress  is  assumed  to  be  r,  then  from 
equation  (4) 


5 


NACA  TN  2692 


1  + 


n+1 

1 

nv5n 


=  r 


(5) 


and  the  relative  thickness  of  the  laminar  sublayer  is 


’L  _ 


n(r  -  l)  v~| 

I 


n 

n+l 


u 


(6) 


The  velocity  in  the  boundary  layer  at  this  point  may  be  computed  from 
equation  .  (2)  and  is  found  to  be 


UL 


u. 


n(r  -  l)  V 


k2  u06 


1 

n+1 


(T) 


since 


In  the  region  0  ^  y  ^  the  flow  is  principally  laminar  and, 

=  p  ~  =  0  at  the  wall  as  may  be  seen  from  equation  (l), 

Sy  dy2 

the  velocity  profile  is  essentially  a  straight  line  in  this  region. 

The  boundary  layer  is  now  assumed,  for  the  purpose. of  determining  the 
wall  stress,  to  be  sharply  divided  into  a  turbulent  region  having  a 
power-law  velocity  profile  and  a  laminar  region  having  a  straight-line 
velocity  profile  as  shown  in  figure  1  by  the  solid  lines.  The  real 
state  of  affairs  is  indicated  approximately  by  the  dashed  line.  The 
stress  in  the  laminar  layer  just  below  the  assumed  intersection,  the 


point  where  6  =  Sp  and  u  =  uT ,  is  now  given  by  (\ 1  — ) 

V  dy/L 


which  will 


be  taken  as 


HUL/5L. 


Since 


Sy 


0  at  the  wall  or  the  stress  in  the 


laminar  sublayer  is  approximately  constant. 


Tw  5. 


1 — 1 

1 

u 

LA 

6 

L  k2 

u0S 


1-n 

n+1 


(8) 
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The  skin-friction  coefficient  is 


41 


1-n  2 


For  n  =  7,  equation  (9)  becomes 


cf  = 


2 


k2 

7(r  -  1) 


1 


=  Constant  R 


(10) 


Equation  (10)  is  one  form  of  the  skin-friction  law  that  has  been 
shown  to  be  very  close  to  the  true  state  of  affairs  for  the  turbulent 
boundary  layer.  The  value  of  the  constant  found  experimentally  is  0.045. 

2 

(See  reference  1,  p.  14-7. )  The  value  of  the  ratio  — ^ — .  which  will 
yield  this  result,  is  0.0444. 


r  -  1 


The  assumption  of  a  definitive  ratio  of  .stresses  is  tantamount  to 
assuming  that  the  Reynolds  number  which  is  characteristic  of 

the  laminar  sublayer,  is  a  constant.  This  Reynolds  number  formed  from 
the  values  of  5^  and  UL  from  equations  (6)  and  (7)  is 


_  uL6L  n(r  -  1) 
RL  =  —  =  "  — 


v2 

Putting  n  =  7  and  - 

r  -  1 


0.0444  into  equation  (ll)  gives 


(n) 


rl  =  158 


It  is  interesting  to  consider  the  significance  of  this  number.  It 
may  well  be  close  to  the  critical  Reynolds  number  of  the  straight-line 
velocity  profile  in  the  laminar  sublayer.  Below  such  a  Reynolds  number  * 

all  disturbances  would  be  damped;  whereas  above  such  a  Reynolds  number 
disturbances  would  feed  on  energy  from  the  steady  motion.  This  condition 
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would  define  the  laminar  sublayer  thickness  and  also  a.  mechanism  by 
which  the  turbulent  boundary  layer  might  receive  energy  in  order  to 
continue  its  existence. 

An  estimate  of  this  stability  criterion  can  be  made  by  using  the 
results  of  von  Karman  (reference  9)*  The  size  of  roughness  height  h 
required  just  to  start  to  change  the  skin  friction  was  found  to  be  given 
approximately  by 


also 


h  ^  1 
SL  ~  4 


If 


'T  ^ 

1  w  ~ 


^UL 

8l  ' 


ULSL 

. 

v  ~ 


lkk 


The  agreement  between  this  value  and  the  value  obtained  by  the  present 
analysis  is  good. 


An  idea  of  the  accuracy  of  the  boundary- layer  model  assumed  may 
be  had  by  referring  to  figure  2,  which  is  a  replot  of  the  data  pre¬ 
sented  in  figure  5  of  reference  10.  The  velocity  profile  shown  is  that 
in  a  fully  developed  turbulent  channel  flow  where  the  Reynolds  number 
based  on  the  half  width  of  the  channel  is  12,200.  The  sharp  transition 
from  laminar  profile  to  turbulent  profile  is  quite  evident.  If  the 

value  for  — found  for  the  condition  of  n  =  7  is  substituted 
r  -  1 

into  equation  (6),  there  results 


7 


(12) 


For  a  value  of  R^  =  12,200,  equation  (12)  yields  a  value  of  —  =  0.022 

which  agrees  well  with  the  experimental  value  of  0.025  determined  from 
the  intersection  of  the  straight  lines  in  figure  2. 
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EXTENSION  TO  COMPRESSIBLE  FLOWS 


The  extension  of  the  foregoing  derivation  to  compressible  flows 
is  straightforward.  Equation  (4)  for  the  ratio  of  the  total  stress  to 
the  laminar  stress  still  holds;  however,  since  the  value  of  v  is  not 
constant,  it  must  take  on  its  local  value  dependent  on  y.  The  ratio 
of  stresses  at  the  edge  of  the  laminar  layer  is  therefore 


so  that 


and 


T 

Tlam 


r 


n+1 

u0k25L  n 

1  +  -2 - h - 

.  1 
nv^5n 


n(r  -  l) 
k2 


VL 

uq5 


n 

n+1 


uo 


n(r  -  l) 
k2 


vLl  n+1 


(13) 


(14) 


When  the  stress  at  the  wall  is  evaluated  it  must  be  kept  in  mind 
that  the  value  of  the  viscosity  is  not  constant  through  the  laminar 
sublayer.  The  velocity-profile  shape  assumed  is  shown  in  figure  3* 

The  stress  in  the  laminar  sublayer  just  below  the  intersection  is  again 

I  ,  which  will  be  taken  as  Although 

the  viscosity  may  vary  as  the  wall  is  approached,  the  stress  must  remain 
constant  as  before;  therefore  the  resulting  wall  stress  is 


assumed  to  be  given  by 


/  du 
V 


_  _  ^LUL 

_  ^Luo[i 

i(r  -  1)  VL 

n+1 

w  "  sL 

5 

k^  u06 

The  skin- friction  coefficient  becomes 


1-n  2 


n(r  -  1) 

n+i/vL  ' 

\n+1  PL 

k2  „ 

W5, 

Po 

(15) 


(16) 
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When  equation  (l6)  is  written  in  the  form 


(17) 


it  may  be  compared  with  the  value  for  incompressible  flow  given  by 


If  the  value  of  k  is  not  changed  appreciably  in  going  to  fairly 
high  Mach  numbers,  for  example,  0  M  ^  4  or  5*  then  for  the  same  value 
of  and  n  in  each  case. 


cf 

CfM=0 


2 


fL 

Po 


(18) 


The  magnitude  of  the  reduction  in  skin  friction  with  Mach  number 
as  given  by  e quat ion • ( 17 )  can  be  found  in  the  following  manner.  Since 
the  pressure  is  constant  throughout  the  boundary  layer 


PL  = 

Po  Tl 

and  when  it  is  assumed  that 


(19) 


n 

Po 


equation  (l8)  becomes 


(20) 


(21) 
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Since  the  discussion  is  for  zero  heat  transfer  and  since  the  assump¬ 
tion  of  constant  total  energy  in  the  turbulent  layer  is  not  liable  to 
introduce  too  much  error,  the  temperature  ratio  in  equation  (2l)  may  be 
written 


(22) 


Thus, 


cf 


c 


fM=0 


1  + 


7  -  1 

2 


M‘ 


_  l+2m-n 


(23) 


For  air  where  7  =  1.4,  m  =  0.76,  and 


n  =  7)  equation  (23)  becomes 


Equation  (24)  shows  that  the  ratio  of  skin-friction  coefficient  for 
the  compressible  case  to  the  incompressible  skin-friction  coefficient 
depends  principally  on  Mach  number  and  only  to  a  very  limited  extent  on 


Reynolds  number  since  ul/u0  is  not 
number. 

The  value  of  f— ^  to  be  put 

Vuo/ 


equation  (13)  as 

follows: 

1 

= 

n(r  -  1)  vL 

n+1 

uo 

k2  uo5 

particularly  sensitive  to  Reynolds 


into  equation  (24)  is  found  from 


1  1 
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2  2(hh-1) 

n(r  -  1)  Vo1n+1/TL\  n+l 


Thus, 


n(r  - 

i^2 


vo 

n+1 

u0S 

1  +  1 


Vuo 


-»  2(m+l) 
2]  n+1 


where 


vo" 

m+1 

u0^ 

Equation  (25)  may  he  solved  for  (— \  for  any  Mach  number  and 

Wo  / 

Reynolds  number  R§.  These  results  may  then  be  inserted  into  equa 
tion  (23)  for.  Cfjcf  . 


»l\2 


Figure  4  is  a  plot  of  j  and  C^/C%  0  aSa^-ns't  Macil  number 

obtained  from  equations  (24)  and  (25)  for  several  typical  values  of  R5 

for  a  one-seventh-power  velocity  profile.  The  insensitivity  of  skin- 
friction-coefficient  ratio  to  Reynolds  number  is  quite  evident. 

It  is  always  interesting  for  the  sake  of* argument  to  see  the  effect 
on  the  skin-friction  ratio  as  the  Mach  number  approaches  infinity, 
although  in  this  case  such  an  extension  is  not  justifiable.  As  the  Mach 
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/ul\2  , 

number  approaches  infinity,  the  term  )  approaches  unity  (see 

equation  (25) ),  hut  the  expression  inside  the  brace  in  equation  (24) 
approaches  l/A,  so  that 


n-l-2m 


lim 
M— >00 


/  cf\_ 

n(r  -  1)  vG 

\CfM=0/ 

k2  uo6 

(n+l) (m+l) 


COMPARISON  WITH  EXPERIMENT  AND  OTHER  METHODS 


In  figure  5  the  results  of  this  analysis  are  compared  with  the 
experimental  results  presented  by  Wilson  (reference  7)  and  with  the 
experimental  and  theoretical  results  of  Rubesin,  Maydew,  and  Varga 
(reference  6)  as  well  as  with  the  original  suggestion  by  von  Karman 
(reference  ll).  The  result  obtained  theoretically  by  Van  Driest  (refer¬ 
ence  8)  for  a  Reynolds  number  Rx  equal  to  7  X  10^  is  also  compared. 
This  result  may  be  compared  since,  the  values  of  R^  used  in  figure  5 

cover  a  range  sufficient  to  encompass  the  value  of  5  that  might  exist 
at  x  in  the  analysis  of  Van  Driest.  The  agreement  between  this  method 
and  the  experimental  results  to  date  which  are  described  very  well  by 
the  curves  of  Rubesin  and  Wilson  is  good.  The  curve  obtained  by  the 
method  suggested  by  Monaghan  (reference  4)  almost  coincides  with  the 
experimental  curve  of  Wilson  and  therefore  is  in  good  agreement  with 
that  of  the  present  method. 


CONCLUSIONS 


A  derivation  of  the  form  of  the  incompressible  turbulent  skin- 
friction  law  for  an  insulated  flat  plate  is  made  in  such  a  way  that  it 
may  be  extended  to  compressible  flows.  The  following  conclusions  may 
be  drawn: 

(l)  The  incompressible  argument  demonstrates  that  a  Reynolds  number 
characteristic  of  the  laminar  sublayer  is  a  constant  approxi¬ 

mately  equal  to  153.  This  value  may  be  close  to  the  critical  Reynolds 
number  of  the  straight-line  velocity  profile  in  the  laminar  sublayer 
below  which  all  disturbances  are  damped  and  above  which  disturbances 
may  feed  on  the  energy  of  the  steady  motion. 
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(2)  An  expression  for  the  ratio  of  compressible  skin-friction 
coefficient  to  incompressible  skin-friction  coefficient  at  the  same* 
value  of  Reynolds  number  based  on  the  boundary- layer  thickness  on  an 
insulated  flat  plate  is  derived.  '  The  agreement  between  the  compres¬ 
sible  and  incompressible  skin-friction  coefficients  found  experimentally 
and  those  predicted  by  the  present  method  is  good. 

(3)  The  ratio  of  compressible  to  incompressible  skin  friction  is 
shown  to  be  rather  insensitive  to  Reynolds  number  based  on  the  boundary- 
layer  thickness  and  to  be  principally  a  function  of  Mach  number. 
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Present  method 


Figure  5.-  Comparison  of  experiment  with  present  method 
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